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Abstract 



On 
O 

Os 

We use Heavy Quark Effective Theory (HQET) techniques to parametrize certain 

Æ ■ 

non-perturbative effects related to quantum fluctuations that put both heavy quark and 



^ ■ antiquark in quarkonium almost on shell. The large off-shell momentum contributions are 
calculated using Coulomb type states. The almost on-shell momentum contributions are 
evaluated using an effective 'chiral' lagrangian which incorporates the relevant symmetries 
of the HQET for quarks and antiquarks. The cut-off dependence of both contributions 
matches perfectly. The decay constants and the matrix elements of bilinear currents at zero 
recoil are calculated. The new non-perturbative contributions from the on-shell region are 
parametrized by a single constant. They turn out to be 0(a 2 /ÅQCDan), &n being the Bohr 
radius and a the strong coupling constant, times the non-perturbative contribution coming 
from the multipole expansion (gluon condensate). We discuss the physical applications to 
T, J/ty and B c systems. 



1. Introduction 

The so-called Heavy Quark Effective Theory (HQET) [1-5] has become a standard tool 
to study the properties of hadrons containing a single heavy quark (see [6] for reviews). 
The hadron momentum is essentially the momentum of the heavy quark which may then 
be considered almost on-shell. The dynamics becomes independent of the spin and the 
mass of the heavy quark giving rise to the so-called Isgur-Wise symmetries [1,2]. The 
relevant modes are momentum fluctuations of the order of Aqcd which are described by 
the HQET [3-5]. One cannot actually carry out reliable perturbative calculations at that 
scale, but one can certainly use the Isgur-Wise symmetries to obtain relations between 
physical observables. 

For hadrons containing two heavy quarks or more the HQET is not believed to be a 
suitable approximation. The reason being that a system of two heavy quarks is mainly gov- 
erned by the perturbative Coulomb-type interaction. The relevant modes are momentum 
fluctuations of the order of the invers Bohr radius, which is flavor dependent, and not of 
the order of Aqcd- Still, if one is interested in subleading non-perturbative contributions 
related to the " on-shellness" of the heavy quarks, the HQET may provide some useful infor- 
mation. Irrespectively of the above, the HQET has already been used in phenomenological 
approaches to two heavy quark systems [7] . 

We shall argue that the leading non-perturbative contributions in the on-shell region 
to the quarkonium decay constants and to the matrix elements of heavy-heavy currents 
between quarkonia states can be described by a suitably modified HQET. The well-known 
non-perturbative contributions in the off-shell region arising from the multipole expansion 
[8,9] are 0(ÅQCDa n /a 2 ), a n being the Bohr radius and a the strong coupling constant, 
times the contributions we find. The key observation is that when the heavy quarks 
are almost on-shell the non-perturbative effects must be important. In that regime the 
multipole expansion breaks down, but it is precisely there where HQET techniques become 
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applicable. 

In ref. [10] it was pointed out that when fields describing both heavy quarks and 
heavy antiquarks with the same velocity are included in the HQET lagrangian, the latter 
has extra symmetries beyond the well known flavor and spin symmetries [1,2]. In ref. 
[11] the extra symmetries were thoroughly analysed (see [12] for related elaborations). It 
was shown that they are spontaneously broken down to the spin and flavor symmetries, 
even if the gluons are switched off. The Goldstone modes turn out to be two particle 
states with the quantum numbers of s-wave quarkonia. Translating these findings into 
phenomenologically useful statements was the original motivation of this work. 

The main hypothesis in what follows is that whenever we have a heavy quark field we 
may split it in two momentum regimes. The momentum regime where the heavy quark 
is almost on shell (small relative three momentum), and the momentum regime where the 
heavy quark is off shell (large relative three momentum). The main observation is that the 
HQET should always be a good approximation for a heavy quark in the almost on-shell 
momentum regime of QCD [10,12], no matter whether the heavy quark is accompained 
by another heavy quark in the hadron or not. What makes a hadron containing a single 
heavy quark qualitatively different from a hadron containing, say, two heavy quarks are 
the large off-shell momentum effects. In the former the large off-shell momentum effects 
are small and can be evaluated order by order in QCD perturbation theory [1,5,13,14]. In 
the latter the large off-shell momentum effects are dominant giving rise to Coulomb-type 
bound states. However, once this is tåken into account there is no a priori reason not to 
use HQET in the almost on-shell momentum regime for systems with two heavy quarks. 
Then the extra symmetries found in [10,11], which naturally involve quarkonium systems, 
should be relevant. 

Suppose we have two quarks Q and Q' which are sufficiently heavy so that the for- 
malism below can be readily applicable. Let us denote by tpQ, tjq, Qq, and Qqi the vector 
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QQ, pseudoscalar QQ, vector QQ' and pseudoscalar QQ' states. Our main results follow. 

(i) The masses do not receive new non-perturbative contribution from the on-shell 
momentum region. Consequently, the leading non-perturbative correction comes from the 
multipole expansion [8,9]. This allows to extract mq in a model independent way from 
m^ Q , and hence fix the parameter Å relating mq with the mass of the Qq systems [6] . 

(ii) The new non-perturbative effects from the on-shell momentum region in the de- 
cay constants /ø Q , f VQ , fq* t and fq , are given in terms of a single non-perturbative 
parameter fu- 

(iii) The new non-perturbative effects from the on-shell momentum region in the 
matrix elements of bilinear heavy quark currents at zero recoil are given in terms of the 
same non-perturbative parameter fjj. In particular, this implies that the semileptonic 

decays (mq > m Q >) 

i>Q, Vq — ► Qq> , Qq' 

Qq' > Qq' — ► rfQ' » VQ' 

at zero recoil are known in terms of fy Q , f VQ , fq* ; and fq Ql - 

We distribute the paper as follows. In sect. 2 we perform so me short distances 
calculations in the kinematical region we are interested in. In sect. 3 we summarize the 
main results of ref. [11] and match the results from sect. 2 with the HQET. In sect. 4 
we construct a hadronic effective lagrangian for on-shell modes in quarkonium. In sect. 
5 we calculate the decay constant. In sect. 6 we calculate the matrix elements of any 
bilinear heavy quark current between quarkonia states. This is relevant for the study of 
semileptonic decays at zero recoil. In sect. 7 we briefly discuss the possible use of our 
formalism for T, B c , B*, J/ty and 1] C physics. Section 8 is devoted to the conclusions. In 
Appendix A we show how to include l/m corrections in the hadronic effective lagrangian 
for the on-shell modes. A few technical details are relegated to Appendix B. 
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2. Short distance contributions in the on-shell momentum regime 



As mentioned in the introduction, what makes a QQ system qualitatively different 
from a Qq system are the short distance contributions. In a Qq system these are well 
understood. They amount to Wilson coefficients in the currents and in the operators of the 
HQET lagrangian, with anomalous dimensions which are computable in the loop expansion 
of QCD. For a QQ system the short distance contributions cannot be accounted for by 
just anomalous dimensions in Wilson coefficients. Indeed, the anomalous dimension of a 
current containing a heavy quark field and a heavy antiquark field with the same velocity 
becomes imaginary and infinite [15]. For large rriQ, the two quarks in a QQ system appear 
to be very close. Due to assymptotic freedom the system can be understood in a first 
approximation as a Coulomb-type bound state. In perturbation theory this is equivalent 
to sum up an infinite set of diagrams (ladder approximation) whose kernel is the tree level 
one gluon exchange (see [16] for a review). 

We shall assume that the dominant short distance contribution to heavy quarkonia 
is the existence of Coulomb-type bound states. Typically we shall be interested in Green 
functions of the kind 

Gripum) ■■= f rf 4 æi d 4 a; 2 e^ æi+ ^^(0|T{Q a rg 6 (0)Q^H^i)QaJ(^)} |0> , (2.1) 
for the range of momentum 

Pl = -m b v -k x , p 2 = -m a v - k 2 , (2.2) 

k\ and k 2 being small. 

Since the quarks are very massive, for the range of momentum (2.2) the leading 
contribution to (2.1) is only given by the following or der ing 

G r (pi,P2)= / d 4 x 1 d 4 x 2 e tpiXl+tp2X2 6(-max(x 1 ,x° 2 )) 

J (2.3) 

x (0|Q a rg b (0)T{gS( æi )QS(a; 2 )}|0). 
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We insert the identity between the current and the fields and we approximate it by the vac- 
uum plus the Coulomb-type states (the states above threshold shall not give contribution 
when we sit in the relevant pole). We treat then the fields as being free. 

/d?P - - 

j^^\s,P n = m ab}n v)(s,P n = m ab , n v\ (2.4) 

The Coulomb state in the center of mass frame (CM) reads 
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-> .1 m\ n f d k ~ /r >. 1 

\S, P n = m a b n V) = , O0- V u / r^W 5 n (fe) — . 

1 ab ' n ' JW c m ah , n J (2tt)3 ab > ny J ^2p°2$ (25) 



xJ2u a (pi)r s vP(p 2 )aUpi)bl( P 2)\0) , 

a ,f3 



where 



p x = m a v + k + — — ttV , p 2 = m b v -k - — — --v , 
1 + v u 1 + v u 

Pi = m a v° + k.v , P2 = mbV° — k.v , 

m ab :=m a + m b , m a6 , n := m ab - E abj7l ,T S = i-y 5 p_,iftp_ , 

v 2 = l , p ± :=^- , e\v = 0. (2.6) 

E a b,m ^ab,n(x) and \& a b,n(k) are the energy, the coordinate space wave function and the 
momentum space wave function of a Coulomb-type state with principal quantum number 
n. v is the bound state 4-vector velocity. a^ a (pi) and b^(p2) are creation operators of 
particles and anti-particles respectively. u a (pi) and v^{p2) are spinors normalized in such 
a way that in the large m limit the following holds 

Y,n a ipi)u a {p l )=p + , '£ i v a (p 1 )v a (p 1 ) = -p-. ( 2 - 7 ) 

a a 

Choosing the momenta as in (2.6) is crucial in order to take into account that the CM 
of the bound state moves with a fix velocity v with respect to the laboratory frame [17]. 
(2.5) has the usual relativistic normalization 

(s,P n = m abyn v\r,P m = rriah^v'} = 2m ab , n v°(27r) 3 5 (3) (m ab , n (v - v'))S nm 5 rs . (2.8) 
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We have to consider the following kind of matrix elements 

{a, m ab , n v\Q a a2 (x 2 )Q b ai ( Xl )\0) = e im ^ v - x (s, m ab , n v\Q a a2 (x 2 - X)Q b ai ( Xl - X)\0) 
_ p im ab , n v.x Kb (f x f d s k ~ -> i( kMx°-x(k +1 ^v)) ( 2 - 9 ) 

m a x x + m b x 2 

X = , x = x\ — x 2 ■ 

m ab 

where it is essential to extract the CM dependence in the fields before using the explicit 
expression (2.5) for the calculation of (2.9). As mention above the states \s,m ab , n v) have 
the explicit expresion (2.5) only in the CM frame [16,17]. Factors of the kind m ab /m ab ,n 
appearing in several expressions above have been approximated to 1 in the rest of the 
paper. Finally, performing the x±, x 2 integral and taking into account that 

-2(P+) a 2 a 3 

(P-) (2.10) 

s 

we obtain 



Gt(pi,p2) = ^*«&,»(o)*a& ) n(o)(p-rp + ) aaai ^ 

n 

1 1 



«2 



(2.11) 



x 



v.k 2 + ^-E ab , n + ie v.h + ^-Ea^n + ie ' 

In the last expression we approximated W ab;n (e\/c) ~ * a 6,n(0) (we neglect 0(( n ^ Q fc ^ ) 2 )). 
In (2.11) there is a sum over an infinite number of poles. Each term in the sum corresponds 
to a Coulomb-type bound state. At the hadronic level we want to describe only one of 
those states. This is achieved by tunning the external momenta to sit on the relevant pole. 
Suppose we are interested in ipQ (ra) state. Then we take 

h = k[ - ^-E abn v , k 2 = k' 2 - ^-E ab n v , (2.12) 

m ab m ab 

so that in the limit k\ — > 0, (i = 1, 2) we obtain 

Gr(pi,P2) =tt* 6) j0)tt o &, n (0)(p_rp + ) a2ai <y ilia 

1 1 (2-13) 



x 



v.k' 2 + ie v.k[ + ie 
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Notice from (2.2) and (2.12) that we must subtract from the momentum of the quark 
(m a — IIh ^E ab n )v in order to get an expression suitable to be reproduced in the HOET. This 
may be interpreted as if integrating out off-shell short distance degrees of freedom produces 
an effective mass for the almost on-shell modes of a heavy quark inside quarkonium. This 
effective mass depends on the precise bound state the quark is in. We are almost on-shell 
when v.k^e 3 .k\ ~ A QC d 0=1,2). 

This restricts the vality of our approximation to the case E ab:Tl ~ V-abO? /n 2 3> A-qcd 
( Hab is the reduced mass), otherwise momentum fluctuations of the order of Aqcd would 
take us from one pole to another. Notice also that for arbitrary large but fix [i ab there is 
always an n where this approximation fails. Therefore we shall always be dealing with a 
finite number of low laying energy levels. 



Consider the four-point function. 

G( Pl ,p 2 ,P3,P4) := j d 4 x 1 d 4 x 2 d 4 x 3 d 4 x 4 e ipiXl+ip2X2+ip3X3+ip ^ 

x (o\T{Q b :i(x 1 )Qznx2)QZ i (^)Q b ::(x,)} io> 



(2.14) 



For the momenta 



Pl = 


-(m b 


E a b,n)V k 1 , 

m ab 


P2 = 


(m a - 


Tfl 

E ab ,n)v + k 2 , 

m ab 


P3 = 


-{m a 


m ab 


P4 = 


(m b - 


E ab ,n)v + k A , 

m ab 



(2.15) 



(k^ — > , i = 1, 4) working in the same way we obtain 

G( Pl ,P2,PS,P4) = (27T)V 4 )(-^ + k' 2 -k' 3 + k' 4 )^ Yl (rn) a2 a 4 (f n ) aiQ3 

r n =i~f 5 p_,ifp- 

(2.16) 

We shall see in the next section that (2.13) and (2.16) can be reproduced (with suitable 
changes) by a HQET for quarks and antiquarks. 
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3. HQET for quarks and antiquarks 



The lagrangian of the HQET for quarks and antiquarks moving at the same velocity 
V/j, (v^v^ = 1) reads [4] 

L v = ihy^v^D^hy = ihyV ■ Dh^ — ih~v ■ Dh~ , (3-1) 

where h v = hy + h~ and hf = ■^-hy. hy contains annihilation operators of quarks with 
small momentum about mv^ and h~ contains creation operators of anti-quarks again with 
small momentum about mv^. is the covariant derivative containing the gluon neid. 

The quark and antiquark sector of (3.1) are independently invariant under the well- 
known spin and flavour symmetry [1,2,4] 

hf -> e l£ ± s ?ht and hf -> h±e~ ie ± s t , (3.2) 

where Sf = i€ijk[$j,$k]0- =t^)/2, with ,j = 1,2,3 being an orthonormal set of space 
like vectors orthogonal to , and 

hf -> e lø± hf and hf -> h±e~ ie± . (3.3) 

e± and 6± are arbitrary real numbers corresponding to the parameters of the transforma- 
tions. 

The lagrangian (3.1) is also invariant under the following set of transformations 

h v e^ e , (3.4) 

h v e^ e , (3.5) 

h v e e *^ , (3.6) 

h v e ié ^ . (3.7) 

The whole set of transformations (3.2)-(3.7) corresponds to a U(4) symmetry for a single 
flavour. For N^f heavy flavours they correspond to a U{ANhf) group. In the latter case 
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hu ^ g hy • hy 



h"u ^ 6 h/y <j h/y 
hu ^ €> ^ hu ^ hu 
h v -> e i£% ^h v ; h v 



h v must be considered a vector in flavour space and the parameters of the transformations 
(3.2)-(3.7) as hermitian matrices in that space. 

When the gluons are switched off it is easy to prove that the U{AN h f) symmetry breaks 
spontaneously down to U(2Nhf) <g> U(2Nf l f) (see [11]). The following currents correspond 
to the broken generators 



a, 6, c... = 1, ...Nhf are flavour indices. They transform according to two four dimensional 
irreducible representations of U (2Nh f)®U (2Nh f). In what follows we are going to assume 
that the situation above is not modified when soft gluons are switched on. The currents 
(3.8) have the quantum numbers of pseudoescalar and vector quarkonium respectively. The 
heavy quark and antiquark fields interact with soft gluons according to the lagrangian (3.1). 
For soft gluons, perturbation theory cannot be realiable applied. However, one can use 
effective lagrangian techniques, which fully exploite the symmetries above, to parametrize 
the non-perturbative contributions in this region. This shall be done in section 4. 

For further purposes let us carry out some leading order perturbative calculations. 
Consider first 



where \x is an ultraviolet symmetric cut-off in three-momentum (see [11] for more details). 
Consider also 



jf ± := h a v i l5 p±h b v and jf± := h a v i^p ± H 



'V 1 



b 



(3.8) 




(3.9) 




(o|T{/irr // ^ + (x 1 )^+r^+(o)^+r'/ir(^2)}|o) 



(3.10) 
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The flavor indices (a , b , c) are not summed up unless otherwise indicated. Colour indeces 
are not explicitly displayed in the colour singlet currents. Otherwise they will be denoted 
by h , ii-, ■■■= 1----N C , N c being the number of colours. We shall drop the subscript v from 
h v and change the superscripts ± into subscripts in the following. 

The strong cut-off dependence of (3.9)- (3. 10) is puzzling. We shall see later on that 
it cancels against suitable short distance contributions. 

As claimed before, it is easy to see that (2.13) is reproduced by the following HQET 
Green function at tree level 

GMX) = f ^i^e-^^-^^^OlTjCr^r/i^O)^^!)^^^)} |0> (3.11) 

with Cr being a Wilson coeficient. 

Cv = *: M (0)* aM (0) . (3.12) 

Analogously, (2.16) is reproduced in the HQET by * 

* One may be tempted to include (3.13) as a perturbation in the HQET lagrangian. 
This is not quite correct. The Green function 

G(k[,k' 2 ,k' 3 ,ki) = J d 4 x 1 d 4 x 2 d 4 x 3 d 4 xie- lk i Xl+ik 2 X2 -< X3+ik * X4 

x (0\T {h^Mlfå (x 2 )KZMh-^M} 1°) 

gives a non-zero contribution in the HQET which does not correspond to (2.14)-(2.16). It 
is (3.13) which gives the leading contribution to (2.14) in the HQET and hence the last 
term in (3.13) must not be included in the lagrangian. This means that unlike in the case 
of heavy-light systems, the short distance effects here cannot always be accounted for by 
only modifications of the currents and the lagrangian, as we may have naively expected. 
We have to content ourselves by identifying for a given Green function, the Green function 
in the HQET which gives the same result. 
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G(k[, k' 2 , k' 3 , k' 4 ) = J rf 4 xirf 4 a:2rf 4 a;3rf 4 X4e- jfc i Xl+ ^^-^3-3+ l ^x 4 
x (0^1^^)^ (x 2 )h a £ a (z 3 )h b _^(zt) 
xijd 4 y (-^*:^(0)* a6 , n (0)) ^ 

(3.13) 

4. Effective hadronic lagrangian for the on-shell contributions of s-wave 
quarkonia 

We have seen that for the on-shell kinematical regime certain correlators can be re- 
produced in the HQET. We shall see in the sect. 5 and 6 that the contributions from this 
region to the decay constants and matrix elements reduce to the evaluation of heavy quark- 
antiquark currents in the HQET. For the range of momentum we are interested in these 
Green functions cannot reliable be evaluated in perturbation theory. We shall use in this 
section effective lagrangian techniques, very similar to those used in Chiral perturbation 
theory, to parametrize the nonperturbative contribution. 

There are well-known rules [18] (see also [19]) to construct phenomenological la- 
grangians for Goldstone bosons associated to the symmetry breaking of a group G down to 
a subgroup H for relativistic theories. These rules need two slight modifications to become 
applicable to our case: 

(i) The HQET is formally relativistic only after assigning transformation properties 
to the fix velocity v^. We must take into account that the velocity as well as the 
can also be used to build up relativistic invariant terms. 

(ii) The HQET is not only globally U(4N hf ) invariant, but locally U(4N hf ) gauge 
invariant under transformations which only depend on the components x 1 := x^é 1 ^. We 
shall also require the phenomenological lagrangian to be local gauge invariant under the 
corresponding transformations. 
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With the above modifications (i) and (ii) we shall apply the rules [18] to the case 
G= U(4N hf ), H= U(2N hf ) <g> U(2N hf ). Let us first associate to the currents (3.8) fields 
in the phenomenological lagrangian which have the same transformation properties under 
H 

H ab -> h a i<y 5P+ h b , H abl -> h a ifp+h b , 

(4-1) 

H ba * ^h b i-f 5 p-h a , H bal ^-h b ifp-h a . 
We build up the following object 

H = i^p-H - i^p-W + i<y 5P+ H^ + ifcp+H lt , (4.2) 

H := 7°H t 7° = H, 

where we use matrix notation for H ab and H ab% . H transforms under the unbroken sub- 
group as follows 

H -> hm' 1 , h e U(2N hf ) <g> U(2N hf ) . (4.3) 

We assign non-linear transformations under the full group U (AN^f) in the standard manner 
[18] 

g(e)e H =: e H '/i(H, 9) , 

geU(W hf ) , heU(2N hf )®U(2N hf ) , e H e U(4N hf )/U(2N hf ) <g> tf(2W fc/ ) , 

(4.4) 

where H' is the transformed neid. Then 

e H — > e H ' = øe 11 /»" 1 = te H S , (4.5) 
where <7 = 7 (7^7°. The following property holds 

^e H =e" H ^, (4.6) 

which implies that 

S := e 2I V = ^e" 2H , S 2 = l , e 2n f ^ ge^fg- 1 . (4.7) 
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Because of the local gauge symmetry we can only build the following connection and 
covariant tensor 

V := \ (e~ n v.de n + e n v.de~ n ) , V — ► hVh' 1 + hv.dh' 1 , fV = Vf, 

A ■= i ( e - K v.de n - e u v.de~^) , A — > hAh' 1 , fA = -Aj> , (4.8) 

v.dS = e R Ae K i>. 

Notice that any derivative with respect to v acting on functions of x % which 

are not scalars will not be covariant under the local transformations. 

The u(4Nhf) algebra and the HQET lagrangian are invariant under the following 
discrete symmetry 

e^-e; , Vp^-v», (4.9) 

which is reminiscent of charge conjugation. They are also invariant under the 50(3) 
transformations — > and, of course, under Lorentz transformations if we assign 

V/j, — > k^v v , — > A^ u e l u . All these symmetries should also be implemented in the 
effective lagrangian. 

We can start at this point the construction of the effective lagrangian, order by order 
in derivatives, using the objects dermed above. At first order it turns out that there is no 
invariant term. Still there is a term which is invariant up to a total derivative. It reads 

Tr(j)V) ~ -Atx{H ] v.dH + H^v.dH 1 ) + ... , 

Tr(i)V) — > TrtøV) + Tr^hv.dh' 1 ) . (4.10) 

Tr means trace over flavour and Dirac indices whereas tr means trace over flavour in- 
dices only. We keep tr for trace over Dirac indices only. It is not difficult to prove that 
Tr^hv.dh -1 ) is indeed a total derivative. This is analogous to the case of the Heisenberg 
ferromagnet where the leading order term in the effective lagrangian for the Goldstone 
mode is also invariant up to a total derivative [20] . Then at leading order the long distance 
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properties of heavy quarkonia are governed by a single constant. At next to leading order 
we have the term 

Tr(AA) ~ -Atx{v.dH ] v.dH + v.dH^v.dW) + ... . (4.11) 

Terms containing x 1 derivatives start appearing at sixth order. Notice that there is no 
vertex involving an odd number of fields. This holds at any order in derivatives and it is a 
consequence of the separate conservation of the number of heavy quarks and antiquarks. 
For convenience we normalize the effective lagrangian as follows 

f 2 + 

-iOL T r{j)V) = ztr(nV<9n + ir Varp) + ... , 

n ir 
h= t - , h* = — . (4.12) 

JH JH 

ffj is a dimension 3 parameter of the order of A.q Cd . The effective lagrangian built above 
makes sense by itself as a toy model. If we ignore the matching with high energies we 
can withdraw some consequences out of the lowest order lagrangian. These and the l/m 
corrections to this toy model are worked out in the appendix B. 

Let us next discuss how to represent quark currents in the effective lagrangian. Con- 

sider 

jf = h a Yh b . (4.13) 

Let us introduce a source af for each of these currents and write all possible currents up 
in the lagrangian 

L v = ihfv^D^h + hfah , 



a:=J2 a r b ^ T - ( 4 - 14 ) 
r 

L is now locally invariant under U(4Nhf) if we assign to a the transformation property 

a — > gag' 1 + giv.dg' 1 . (4.15) 
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At the hadronic level we may also require local gauge invariance upon the introduction of 
a. This is easily achieved by changing v.d into v.d — ia in the definition of V in (4.8). We 
obtain 



e'2 



L = -i J -jL \Tr{j)V) - iTr(aS)} . (4.16) 



Then we may identify 



f 2 

h a Th b — ► -^fTr (TT ab e m ) , (4.17) 



where T ab is the zero matrix in flavor space except for a 1 in row a column b. It is interesting 
to observe that the U(AN h f) symmetry is so large that any bilinear current of the kind 
(4.13) can be written in terms of a generator of the U(4:Nhf) symmetry. This is the actual 
reason why the identification (4.17) does not involve any extra unknown parameter. It is 
analogous to the case of the vector and axial-vector currents in the Chiral Lagrangian [21]. 



Let us next calculate for further convenience the correlators (3.9) and (3.10) in the 
hadronic effective lagrangian. For (3.9) we have 



Grr>(k) = J d 4 xe- lk - x (0\T{h a _Th b + (0)h b + T'h a _(x)} |0> 



-/ 



d 4 xe- lk - x (0\T { -&Tr(p-T P+ T ab e 211( - 0) ) 



-&Tr(p+T'p-T ba e 211 ^) \ |0> 



f 2 

J H ' 



d 4 xe- lk - x (0\T <j -^-Tr(p.Tp + T ab 2R(0)) 



-^i-Tr(p + T'p_T ba 2R(x)) \ 0) 



(4.18) 
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For (3.10) we have 

Gitt» (&£, fea) = / d i x 1 d l x 2 e ik ^ Xl - ik ^ X3 (0|T {h a _T''h b + (x 1 )h\Th c + (0)h c + T'h a _(x 2 )} |0) 

/I 



/ 



x 



d 4 xid 4 x 2 e i *i a,1 - iAi a iB2 (0|T| [-^Tr(p_r> + r 6 e 2H(ll) ) 
_^ Tr(p+rp+T 6c e 2H(o) ) j J_^Tr(p + rV-T ca e 2H ( X2 )) J> 1 0> 

fl, 



J d 4 x 1 d 4 x 2 e< Xl - ik '^(0\T^ ^- J -^Tr(p_T"p + T ab 2R( Xl )) 



x 



^-Tr(p + rp + r te 2H 2 (0))] [-^Tr(p + r'p-r-2H(x 2 )) 



4 tr(p _r-'p + r P+ r-) -^-^- 



(4.19) 



Notice at this point that we may obtain (3.9) and (3.10) from (4.18) and (4.19) by taking 
f]j/2 — > N c fi s /6n 2 . Hence /jj at the hadronic level plays the ro le of the cut-off n at quark 
level. Observe also that the dependence on the T-matrices in (4.18)-(4.19) is explicit. All 
decay constants and matrix elements of bilinear currents are given in terms of the only 
non-perturbative parameter This is a direct consequence of the U(ANhf) symmetry 
being spontaneously broken down to U(2N h f) ® U(2N h f). 

5. Example: the decay constant, f-f 

5.1. Separating and evaluating off-shell and on-shell contributions 

Consider the curr ent- cur rent correlator 

G T {p) := f d 4 xe ipx (0\T {Q a TQ h (Q)Q b TQ a (x)} |0> , (5.1) 
p = -m abjn v-k , k^O. 

We separate 

Q a TQ b = (Q a TQ b ) on + (Q a TQ b ) off . (5.2) 

Where (Q a TQ b ) on and (Q a TQ b ) ff means that both heavy quark fields in the current have 
momenta almost on-shell and off-shell respectively. Our goal is to obtain a representation 



16 



in terms of the HQET of any Green function containing an (Q a TQ b ) on . In order to enforce 
"on-shellness" it is convenient to make the substitution 

j d 4 x(Q b TQ a (x)) on e* px — j d A x 1 Q b ^ i (x 1 )e^ J d 4 x 2 Q a a \ 2 (x 2 )e^ (f ) aia J ili2 , 

(5-3) 

Pl = -(m a -E abjn )v - k[ , 

m ab 

P2 = -(m b - —£„(,,„)!) - k' 2 , 

m ab (5.4) 

k — — k-^ ~\~ k<2 , 

and see whether the new Green function admits a representation in terms of the HQET. 
This is nothing but the calculations carried out above. Then we undo (5.3) by putting the 
fields depending on x\ and x 2 in the HQET at the same point x. We have (from (2.1), 
(2.13) and (3.11)) 

f d 4 xe^(0\T{(Q a FQ b (0)) off (Q b TQ a (x)) on } |0> 
' ^ (5-5) 
d 4 xe- lkx (0\T{C r h a _Fh b + (0)h b + rh a _(x)} |0) . 



Analogously, using (2.14), (2.16) and (3.13) we have 
J d 4 xe ipx (0\T{(Q a TQ b (0)) on (Q b TQ a (x)) on } |0> 

Xijd^y (-2^*:6,n(0)*a6,n(0)) £ ^«^(^iv.^f „fc 6 (y)) ||0> . 

(5.6) 

The contribution involving only off shell quarks has the familiar form 

f d 4 xe^(0\T{(Q a TQ b (0)) off (Q b TQ a (x)) off } |0> 
J 1 (5-7) 

= -tN c tr^P+^P-)\^ab,nm 2 —r—- ■ 

v.k + te 
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The expressions (5.5) and (5.6) correspond to corrections 0{K^- }CD a\ b n ) and 0{Ps^ CD a & ab n ) 
respectively to the leading result (5.7), a a b^ n ~ n/(a/i a b) is the Bohr radius. Since we are 
only interested in the leading non-perturbative corrections we shall neglect (5.6) in the 
following. Let us only remark that the hadronization of the four quark operator in (5.6) 
introduces new parameters. This is because it is not a generator of the U(4Nhf) symmetry 
as the currents of the kind (4.17) are. 

The r.h.s. of (5.5) can be hadronized and calculated using the effective lagrangian 
discussed in section 4. From (4.18) we obtain 

d 4 xe^(0\T{(Q a TQ b (0)) off (Q b TQ a (x)) on }\0) 



'4r(p-Tp + T)*l bin (0)* ab , n (0)f H - 



, (5-8) 



Notice that the result is spin independent and the flavor dependence resides only in the 
wave function, which is known. We finally obtain 

\u Qin) \ 2 = 4m ab , n (N c \^ ab , n (o)\ 2 + l - (*: 6jn (o)* a6 , n (o) + K b , n m^)) fl) , 



(5.9) 

\fipQ(n) I 
1T^ab,n 

Notice that the non-perturbative correction we find to the decay constant is 0{A!q CD a? ah n ) 
and hence presumably more important that the correction arising from the multipole ex- 
pansion which is 0((Aqq DO-ab,n) 4 1 '« 2 ) [8,9] (we count the quark condensate as 0(Aq CD )). 

5.2 Cut-ofF independence 

Let us next discuss the important issue of the cut-off independence. Even though 
we have not written it down explicitely, the introduction of a cut-off to separete almost 
on-shell momenta from off-shell momenta is necessary. Of course, the final results must not 
depend on the particular value of the cut-off. At the short distance end of the calculation, 
the cut-off must exclude momenta which are almost on-shell. This is easily achieved by 
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cutting off small momenta from the wave function 

*«m(0) = / (Ør*«6,»$ — ^ (03*«m(^) == , (5.10) 

where is a symmetric IR cut-off in three momentum. The wave functions in (5.9) must 
be understood as the cut-off wave functions (5.10). On the HQET side the cut-off must 
be ultraviolet. It has already been displayed in the leading order perturbative calculation 
at quark level in section 3. In particular, from (3.9) we obtain 

f d 4 xe^(0\T{(Q a rQ b (0)) off (Q b rQ a (x)) on } |0> 

_■ 3l ( 5 - n ) 

= ^N c tr(p-r P+ r)K b , n m ab , n (o)( p 1 



This strong cut-off dependence, however, is totally compensated by (5.10). Indeed, 
once (5.10) is used we have 



*: M (o)* aM (o) + *: M (o)* aM (o) 



^ 2 



2tt 2 

+ O ((pa ab , n ) 2 ) ) , 



(5.12) 

^ (^,»(0)*a!„(0)(^)) =*: M (0)*aM(0)^ [l + O ((va ab , n )*)] , 

^ (*2n*(°)*«M(0)(^)) =*: M (0)* aM (0)^ [l + 0((^a aM ) 2 )] . 

Notice that the way in which the cut-off dependence cancels is remarkable. The strong 
cut-off dependence of (5.11) was first found in [11]. It was not clear at all which short 
distance contribution it should cancel against. (5.10) gives the solution to that puzzle. It 
is apparent from (5.8) and (5.11) that fjj in the hadronic theory plays the role of the UV 
cut-off in the HQET at quark level. From (5.12) it is clear that the cut-off \i must be 
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much smaller than the invers Bohr radius. Therefore our formalism becomes exact in the 
following situation 

m a , m b > l/a a b, n > A* > ^QCD > k' , 

^ 2 (l/a ab , n ) >>k , 
n l 

Furthermore, we have to assume that (i can be tåken large enough so that we may enter 
the asymptotic freedom regime from the HQET side. Otherwise the matching we have 
carried out at tree level would not make much sense. 

From the discussion above it should also be clear that (5.9) can be written in a cut-off 
independent way at 0((fia a b t n) 3 ) by just replacing 

fu — Ph ■= fl - , (5-14) 

where f H need not be positive. 

5.3. Physical state normalization 

There is still a subtle point which makes eq. (5.9) with the replacement (5.14) not 
quite correct. It has to do with the normalization of physical states. It will be clear later 
on (see eq. (6.14) below) that the states we obtain by this procedure do not have the 
standard relativistic normalization as they are supposed to. When we evaluate the Green 
function (5.1) we insert resolutions of the identity which are approximated by Coulomb- 
type states. This is all right. However the low momentum tale of these states is cut-off 
and substituted by a quantity evaluated using the effective hadronic theory. After doing 
so there is no guarantee that the resolution of the identity we introduced is still properly 
normalized. This can be fixed up by changing 

?/ wk*J n)(nl wk^ n)(4p)N ^' M (5 ' 15) 

where \n){n\^ symbolises the cut-off Coulomb states whose low energy tale is evaluated 
in the hadronic effective theory. We present a heuristic calculation of N n (/j,, /#). 
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We start from the Coulomb-type bound state (2.5) and separate high and low relative 
momentum according to 



|r n , P n — iTlab,nV) — |r n , P n — T^ab,nV) ^ ~\~ |Tn5 Pfi — ^ab,nV) ^ 



(5.16) 



The high momentum part of the physical state can be well approximated by the Coulomb- 
type contribution so we may leave it as it stands. However, the low momentum part receives 
non-perturbative corrections, which we evaluate using the effective hadronic lagrangian. 

We proceed as follows. Since a a b,nfJ> < 1, we can approximate the low momentum 
region by 



|J- ri-, 1 n — llb ab,n iJ l — 



^ab,n(0) 



m 



nsr c ^2m a v°2m b v° m ah , 
• fc< " d 3 k 



x 



/ K ^V d 6 k «-^ 
(2^3E^(^) r ^(^)4(Pl)^(P2)|0), 



(5.17) 



Observe now that (5.17) is nothing but the integral of a local HQET current. 

3 

|r n ,P n =m ah , n v) k <» ~ -^=^^^(0) — ^- / d 3 xe-^h a T n h b (x)\0), 



(5.18) 



where k — > and only low momenta are allowed. 

At this point, we can hadronize the current (see (4.17)) and calculate the low momen- 
tum contribution to N n (fi, fn) 

^(Tn? P-n = r Wlab,nV\^ "n? P n = iTlab,nV ) ^ 



= 2m a6 , n t; (27r) 3 5( 3 )(m a6 , n (t;-^))^l^a6,n(0)| 2 . 
Then, putting together high and low momentum contributions, we have 

(Tn? Pn T^ab,n^\^ni P n ^ab^-nV ) 



(5.19) 



= 2m ab;n v°(2n) 3 5^(m ab , n (v-v')) 
= 2m ab , n v°(2n) 3 5^(m ab , n (v-v')) 



d 3 k ,~ f 2 



fc>At 



(2nY 



^abAk)\ 2 + j^\^ab,n(0)\'- 



(5.20) 



Ph 



1 + ^1^(0)1' 
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Where ffj is defined in (5.14) and notice that the result is cut-off independent. 
Finally, the normalization factor reads 

N n (n,f H ) 



l+ifø|*«M(0)| : 

N n (n, fn) can also be obtained from requiring 



(5.21) 



(r n ,P n =m ab , n v\ J d 3 xQ b -f°Q b (x)\T n ,P^ = m ab , n v') = 2m ab , n v (2n) 3 5^ (m ab , n (v - v 1 )) 

(5.22) 

as we shall see later on. Once we have tåken into account the correct normalization (5.9) 
reads 



l/^ Q (n)| 2 — 4m a6jr 



N c \* ab , n (o)\ 2 + \ (*: 6jn (o)* a6 , n (o) + *: M (o)* aM (o)) ji 



|*a6,n(0)| 2 |*aM(0)| 2 A 



(5.23) 



We shall relegate to section 7 the discussion on the applicability of the limit (5.13) 
and formula (5.23) to physical situations. 

6. Matrix elements at zero recoil 

We are interested in Green functions of the kind 

GrrT»(pi,P2) = j d 4 x 1 d 4 x 2 e ipiXl+ip ^(0\T {Q a T"Q b (x 1 )Q b rQ c (0)Q c r'Q a (x 2 )} \0) ■ 

(6.1) 

For the momentum range 

Pi = m abjn v + k[ , p 2 = -m aCjrn v - k' 2 . 



ki , k 2 ► 



(6.2) 
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We separate each current in almost on-shell momenta and off-shell momenta according to 
(5.2). The leading contribution is given by the term 

GrrT"(pi,P2) 

= J dWx 2 e^ 1+ ^ 2 (0|T{(Q a rV(^^ |0) 

/d s k ~ - ~ -> 1 1 

(6.3) 



and the next to leading contribution by the term 

Gttt"(pi,P2) ■= G^ r „(p 1 ,p 2 ) + G^ r „(p 1 ,p 2 ) + G£r,r»(pi,P2) » ( 6 - 4 ) 
G^} r , / (p 1 ,p 2 ) 

= I d^xie^+^iOlT {(Q a T' , Q\x 1 )) on (Q b TQ c (0)) off (Q c T'Q a (x 2 )) off }\0) 1 

(6-5) 

CrrT" (Pi»P2) 



(6.6) 

G rr'r» G>i»P2) 

= y d 4 Xl d 4 x 2 e ip ^ +ip ^(o\T {(Q a r // Q 6 (^i))o//(Q 6 rg c (0)) o// (QT / g a ( æ2 )) on }|0). 

(6.7) 

The calculation of (6.5) and (6.7) is analogous to the ones carried out in section 2. We 
obtain 

G^ r „(pi,p2) = / ^i^^e^^-^^zC^OlTj^r^^O^rp+r'^ (0)} |0) 

J (6.8) 

x / d Q — — , 

J v.q + ie 

/d 3 k ~ -> ~ 
^J3^,„(fe)*ac,m(fc), 



(6.9) 



p—iqxi 

X / d\ 



v.q + ie 
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/d 6 k ~ - ~ 
— y* abn (k)y ac , m (k). 

Notice that (6.5) and (6.7) can not be written in terms of local Green functions in the 
HQET. One propagator must be kept explicit. 

The calculation of (6.6) is more subtle. We describe it in some detail in the Appendix 
B. We obtain 

G° r n r ^ // (p uP2 ) = J éx l( fx 2 é k '^- ik '^c 2 ^\T {htr'%X(x 1 )h b + Th c + (o)h c + T'h a _(x 2 )}\o), 

(6.10) 

c 2 = vi>: Cim (o)vi> a6 , n (o)^; C)m (o)^ a6 , n (o) . 

This term is the only one in (6.4) which remains in the matrix elements (see (6.14) below). 

We calculate (6.8)-(6.10) using the hadronic effective lagrangian (see formulas (4.18) 
and (4.19)). We obtain 

GTrUn,n) - cJ-ftr^T-^r)-^- -^L- , (6.11) 

= cjftr^r-^r)-^ -±- , (6.12) 

GS**) = C 3 f -ftr(p-T"p + Tp + T')-J— - -J— . (6.13) 

2 v.k\ + te v.k 2 + ze 

The matrix element at zero recoil then reads 

(r n ,P n = m ab , n v\Q b TQ c (0)\T m , P m = m aCyn v) 

/ r f 2 \ 

= -y/m ah , n m ac , m tr{t n TT m ) ( J j^K b , n (k)^ ac , m (k) + ^-^ac,m(0)Kb,n( Q ) j > 

(6.14) 

Tn = ilsP- , iftp- for the pseudoscalar and vector particle respectively. The integral in 
(6.14) must be understood with an infrared cut-off fx. From (6.14) it is apparent that our 
physical states are not properly normalized. Indeed, for b = c and V = 7 one should 
obtain (5.22) but one does not. The reason for this has been discussed at the end of sect. 
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5. The solution consist of introducing the normalization factor N n (/j, f H ) defined in (5.21). 
The properly normalized result reads 



{T n ,P n = m abjn v\Q TQ c (0)\T m ,P m = m aC)Tn v} = -^m abjn m aCyrn tr(T n TT m ) 

c2 12 



X 



^K b ,n(^ac,m(k) (l - ^|W(0)| 2 - ^ ac , m '^ 2 



(6.15) 



Notice that the non-perturbative correction depends only on a single parameter j\ 
which may be extracted from the decay constants calculated in section 5. This is a non- 
trivial prediction which turns out to be a direct consequence of the U(4:Nf l f) symmetry 
being spontaneously broken down to U(2N h f) <g> U(2N h f). 



7. Applications 

If the charm and bottom mass were large enough we could apply the results above to 
the physics of T, rj b , B c , B*, J/ty and r\ c . (The top is believed to be too heavy to form 
hadronic bound states and will be ignored). We analyse in this section whether this is 
so or not. In the systems where the formalism actually applies, we are mainly interested 
in estimating the importance of the new non-perturbative contribution rather than in 
obtaining accurate results. The latter is a much harder task which is definitely beyond the 
scope of the present work. 

Let us first focus on bottom. The fact that the almost 'on-shell' momentum excitations 
in heavy quarkonium are Goldstone modes [11] implies that the T and r\h spectrum does not 
received additional non-perturbative contributions. We may then extract the bottom mass 
from the T mass by means of the formulas given in [8,9] which take into account the leading 
order in the multipole expansion. Since we have established a link between quarkonium and 
the HQET we can next use m& to extract Å, the non-perturbative parameter relating the 
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mass of the B- meson to mj. Moreover, taking into account that A is fiavour independent, 
we may next extract the charm mass m c . We summarize the results in the Tab le I. 

In Table I the values we obtain for to& are about a 3% lower than those obtained in 
QCD sum rules [22] but compatible with a recent QCD-based evaluation [23] and with 
the lattice calculation [24]. The values we obtain for Å are a bit lower but otherwise 
compatible with those extracted from QCD sum rules [6]. Our values for m c are again 
about a 6% lower than the typical values in QCD sum rules [22]. We should emphasize 
that our numbers in Table I are model independent. 

We can next extract the non-perturbative parameter ffj from fr (this is done in Table 
II). We use the following formula 



where the 1-loop QCD corrections and the leading correction from the multipole expansion 
[9] * are tåken into account. 

The numbers in Table II are very sensitive to the scale at which a is tåken. Notice 
that we choose a = a(l/cibb,o) in the Bohr radius and binding energy but a = a(mb) 
in the 1-loop perturbative correction included in (7.1). From Table II we see that for 
the actual values of mj and Aqcd = 100, 150MeV the 'on-sheh 1 ' contribution (f H ) does 
not dominate over the condensate, but it is certainly sizeable. For Aqcd = 200MeV all 
corrections are about the same order and for any value of m the 'on-shell' contribution 
dominates over the condensate. 

Observe that the conditions (5.13), in particular ^> /j ^> Kqcd-, m ay be consid- 
ered as reasonable well fulfilled if we take the cut-off /j ~ 700 MeV (see table III below). 

Let us next turn our attention to charm. The charm mass is known not to be heavy 
enough as for the multipole expansion to work [8]. This means that the non-perturbative 
* We use the formula given in ref. [9] which differs from the ones in ref. [8] . 



f r = 2^3771x^66,0(0) 

\ *6fc,o(0)& 
6*66,0(0) 



§66,0(0) f fl 8a(m b ) 



12 3tt 



+ 8.77mg< J B 2 >(^) 6 

2 



(7.1) 
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contribution overwhelms the perturbative one. Therefore any approximation whose leading 
order is a perturbative contribution, like our approach, will not be able to say much about 
charmonium. In particular, for the 'on-shell' contributions the dimculty lies on the second 
last condition in (5.13) being fulfilled. There is little room to accomodate the cut-off \i 
between the invers Bohr radius and Aqqd as should be clear from Table III. We refrain 
ourselves from giving any numbers for charmonium. 

Unfortunately, the situation is not much better for the B c , which has received consid- 
erable attention lately [25-27]. Nonetheless, once we have we shall give some numbers 
in this case in Table IV. 

From Table IV we see that for Aqqd = 100, IbOMeV the contribution of the con- 
densate is too large for the approach to be reliable. For Aqqd = 200MeV we are at the 
boundary of its validity since the 'on-shell' correction is large. We may thus give a rought 
estimate for f Bc only for Aqcd ~ 200MeV, which turns out to be compatible with the 
estimate obtained by QCD sum rules [26], but about a 30% lower than potential model 
estimates [27]. 

From the Table V it follows that the new non-perturbative contribution is not very 
important in the matrix elements between T — B c states. 

The decay constants and matrix elements above receive contributions from corrections 
of several types: 

(i) QCD perturbative corrections to the Coulomb potential (9(a(l/a n )). These have 
been evaluated at one loop level in [28] (see also [23]). 

(ii) Relativistic corrections to the Coulomb potential (9(a(l/a n )) (see also [28,23]). 

(iii) QCD perturbative corrections to the Green functions 0(a(m)). These corrections 
have been tåken into account in (7.1). They correspond to the only QCD corrections in 
heavy-light systems. In our case they are important for the calculation of matrix elements 
at non-zero recoil. 
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(iv) Non-perturbative corrections arising from the multipole expansion in the off-shell 
momentum region 0(AQ CD af l /a 2 (l/a n )) [8,9]. These corrections have also been tåken 
into account in (7.1). 

(v) Finite mass corrections 0(h 2 ^ CD /m) in the hadronic HQET lagrangian. 
8. Conclusions 

We have demostrated that, contrary to the common belief, HQET techniques are also 
useful for the study of systems composed of two heavy quarks. In particular, we have 
identified new nonperturbative contributions to the decay constants and to certain matrix 
elements which are described by a hadronic lagrangian based on the HQET. All these new 
contributions are parametrized at leading order by a single constant fu- This is non-trivial 
and can be traced back to the fact that a U (ANhf) symmetry is spontaneously broken down 
to U{2N hf )®U{2N hf ). 

It is remarkable that strong cut-off dependences coming from a totally different origin 
match perfectly. Indeed, at the off-shell end the cut-off arises from an integral over a 
Coulomb type wave function, whereas at the on-shell end it arises from a Feynman integral. 

We should also stress that we have been able to put in the same context (i.e. the 
HQET) both heavy-heavy and heavy-light systems. This allows for a model independent 
determination of heavy quark masses from quarkonium, which may then be used to extract 
the parameter Å relating the mass of the heavy-light systems to the mass of the heavy 
quark. 

As far as practical applications is concerned, our formalism is suitable for the ground 
state of the T and rji, family. Unfortunately the charm mass is too small for the formalism 
to become applicable in general to J/^> and B c systems. Nevertheless one may stretch it 
in some cases to obtain information on the mass and decay constant of the latter. 
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We have presented a technique which allows to disentangle the on-shell contributions 
from the rest and match them to the HQET. The matching has been carried out at tree 
level. We have already shown in [29] that the matching also goes through at one loop level. 
Nevertheless, a word of caution is needed. It would be desirable to have a more direct and 
systematic derivation of these results from QCD. Progress in this direction is being made 
[30]. 

Let us finally mention that the hadronic HQET lagrangian can easily incorporate 
heavy-light mesons. The formalism can then be extended to the calculation of matrix 
elements between quarkonium and heavy-light systems. The leading non-perturbative 
contributions to those are also given by fjj and another non-perturbative parameter which 
is related to heavy-light decay constants. Non-recoil contributions can also be evaluated 
within the formalism. 

Acknowledgements 

This work has been supported in part by CICYT grant AEN93-0695. A.P. acknowl- 
edges a fellowship from CIRIT. J.S. has benefited from the Fermilab Summer Visitors 
Program. He thanks E. Eichten for illuminating conversations and J.M. Sanchis for send- 
ing his papers in [7] prior to publication. Thanks are also given to J.L. Goity and P. 
Pascual for the critical reading of the manuser ipt. 

Appendix A: a toy model 

Because of the similarity, both in physics and techniques, to the Chiral perturbation 
theory it is interesting to consider a toy model which contains the on-shell contributions 
only. At quark level the model is described by the HQET with quarks and antiquarks 
with the same velocity of sect. 3. At hadronic level is described by the effective hadronic 
lagrangian of sect. 4. 
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Within this model, the interactions between (rjQ,rjQ/), (t)q,^q>), (ipQ,ipQ'), when 
the two particles move roughtly at the same velocity, are described by a single unknown 
constant. This is analogous to the fact that at lowest order in 3-flavour chiral perturbation 
theory the elastic scattering of (7r,7r), (K,K) and (n,K) is also described by a single 
constant. When heavy-light mesons are included in the effective lagrangian the same 
constant describes the elastic scattering of heavy-light mesons with quarkonium. This is 
also analogous to the fact that the local vertex n-ir-N-N at leading order in the Chiral 
Lagrangian is described by the same constant as the (ir, n) elastic scattering. Let us 
mention at this point that when one actually calculates the scattering amplitudes, one 
obtains zero. This has to do with the universality of the leading order effective lagrangians 
for Goldstone modes [18,19,20]. Any theory undergoing a U(4Nhf) spontaneous symmetry 
breaking down to U(2Nhf) ® U(2N} l f) has the same low energy effective lagrangian (4.12) 
provided the rest of the symmetries in the theory are also the same. It was shown in [11], 
that even when the gluons are switched off the spontaneous symmetry breaking occurs in 
the HQET. In that case there is no interaction in the fundamental theory and hence it is not 
surprising that the scattering amplitudes in the effective lagrangian vanish. Universality 
implies that there will be vanishing scattering amplitudes when the gluons are switched 
on as well. 

Within this model one can also treat l/m corrections in a similar way that small 
quarks masses are dealt with in Chiral perturbation theory. At the quark level the leading 
l/m corrections to the HQET are given by a kinetic term 




(Al) 



and a spin breaking term 



Nhf 




Am, 



a 



a=l 




(A.2) 
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The kinetic term (A.l) does not break the global U(2N h f) <g> U(2N h f) symmetry but it 
breaks its local version. In order to construct at the hadronic level terms which break the 
U(4Nf l f) symmetry in the same fashion as (A.l) does, we introduce the w(4A^/ l /)-valued 
sources ø and ai transforming as 

<P — ► g<\>g~ x , 

ai — > gaig~ x + gdig~ x . (A3) 

Then the term 

—dihffydih , 

dih := (Di + ai)h , djij) := Djif — hfai (A4) 

is on one hand invariant under U(ANhf) and on the other reduces to (A.l) upon setting 

/-L. 

/ 2m„ 



ai = , 



1 



2^T b \i>. (A5) 



\ 



At the hadronic level, we must then construct invariant terms linear in ø, which may also 
contain Up to two space derivatives we have 

tr(S<f>) , (A6) 
tr(S(j)diSd % S) , (A7) 
tr(S<j>)tr(diSdiS) , (A8) 

diS := diS + ciiS — Scii . 

We have not written down terms which coincide or vanish upon using (A. 5). 

For the spin breaking term (A. 2) we may introduce a -u(4A^/)-valued source R 1 trans- 
forming as 

R 1 — ► grtg- 1 (A.9) 

so that (A. 2) is substituted by 

hfR l G l h. (A10) 
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We recover (A. 2) upon setting 

R l = 



j)S L . (All) 



4mb 

There are no terms at the hadronic level with the same symmetry transformation properties 
at lower orders in derivatives. The first possible term appears at third order. 

Therefore the leading l/m corrections introduce 3 new parameters. (A. 6) provides a 
mass term 0(A?Q CD /m) and (A.7)-(A.8) give rise to the usual non-relativistic kinetic term. 
The procedure above can easily be extended to any finite order in l/m. 

Appendix B 

We present in this appendix some technical details on the evaluation of the off-shell 
short distance effects carried out in sect. 6. 
Consider the following matrix element 

). (BA) 

Since two different bound states are involved, it is not clear a priori which CM dependence 
one should substract before using (2.5). Nevertheless, translation invariance implies that 

the result of the calculation must fulfill 

(r n , m ab , n v\Q b TQ c (x + a)\T m , m ac , n v) =é^v.a-im ac , m v.a 

x (r n , m abjn v\Q b TQ c ( x) |r m , m aC;n v) 



(B.2) 



We also have 

(r„, m ab , n v\Q b TQ c (x) |r m , m ac , n v) =é m ^v.i-im ac , m v.i 

(B.3) 

x (r n , m ab , n v\Q TQ c (x - £)\T m , m ac , n v) . 
If we assign £ — > £ + a under translations (B.3) fulfils (B.2). If we also require £ to be a 

linear function of x, then necessarily £ = x and the result is well defined. 

(T n ,m abjn v\Q b TQ c ( 

= e ~-- æ -^--- æ (r n ,m aM i;|g 6 rg c (o)|r m ,m aC)n ^ 



_ e !m^„».i-im« !m D.i 



(BA) 

-tr(TnTT m ) I ^K b , n (k)^ ac , m (k)) . 
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Consider next 



(r n ,m abjn v\Q b £(x 3 )Q c £(x A )\r m ,m aCtn v) . (5.5) 



We are in a similar situation as above. However now translation invariance does not fix 
completely the result. Under the same assumptions we obtain 

(^n,m ab)n v\Q b ^(x 3 )Q c ^(x 4 )\T rn ,m aCjrn v) 

_ i(ax 3 + (l-a)x4)(E aCym -E abtn )+im b v.x 3 -im c v.X4 f j_ /-p p \ c 

— c \ N L n)a 4 a ;i "i 3 i 4 (B 6) 



X 



where a is arbitrary and parametrizes the ambiguity. Usually one never runs into cal- 
culations of the kind (B.5) but rather of matrix elements of currents as in (B.l), which 
are not ambiguous. We find expresions like (B.5) in our calculation because we insist in 
enforcing " on-shellness" in certain currents. In our concrete case we have a current with a 
momentum insertion 

Q b TQ c (x)e ip - x , p = (-m b + m c + E ab , n - E ac , m )v . (B.7) 

In order to enforce on-shellness we substitute it by 

%(^3)QS(x4)e*" SB » + ^(r) a , a4 tf i8i4J 

P3 = -(m b E abjn )v-k' 3 , p 4 = (m c -E aCjrn )v + k' 4 , k' 3 , k' 4 -»• , (B.8) 

m ab m ac 

as mentioned in (5.3). However in doing so there is a momentum missmatch (^^■E ab ^ n — 
J ^E ac ^ m )v which should be fixed somehow in order to get back (B.7) in the x 3 = X4 = x 



limit. The most general way of distributing this momentum missmatch between x 3 and x 4 
is by inserting in (B.8) 

e i(p X3+ (l-(3) X4 )(^-E ab , n -^-E ac , m )v _ ^ 

Any (3 is equaly good since we are eventually interested in the limit x 3 = X4 = x. Notice 
that the ambiguity in a in (B.6) is proportional to the ambiguity in j3 in (B.9). Since we 
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can choose (3 at will, we do it in such a way that the dependence in both a and (3 cancels. 
This is how we are able to obtain a representation of (6.6) in terms of the HQET (6.10). 
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Table I. We use Kqcd as an input and take the one loop running coupling constant 
a at the scale of the invers Bohr radius, i.e. a = a(l/a^,o)- F° r the gluon condensate we 
take the fix value < B 2 >= (585Mey) 4 . The error in m& has been tåken from estimations 
of the hiperfine splittings 0{a 2 ), which are also the main source of error in Å. For m c 
the error comes both from Å and the l/m c corrections. The last column gives our model 
independent determination of tob c . 



Table II. We display the relative weight, with its sign, of the 1-loop (a(rrib)), the 
condensate (< B 2 >) and the 'on-shell' {Jh) contribution with respect to the Coulomb 
type contribution (normalized to 1). The last columns display the mass m cr from which 
the 'on-shell' contribution dominates over the condensate and the value of (fn) 1 ^ 3 - 



Table III. We give the cc, bc, bb invers Bohr radius as a function of Aqcd- 



Table IV. We display the analogous to Table II for B c . We have also given our 
predictions for fs c in the last column. 



Table V. We give the relative weight, with its sign, of the 'on-shell' contribution with 
respect to the Coulomb-type contribution (normalized to 1) in the matrix elements (6.15) 
between T — B c states. 
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Table I. 



A QCD (MeV) 


m 6 (MeV) 


A(MeV) 


m c (MeV) 


m Bc (MeV) 


200 


4877±35 


436±35 


1539±70 


6212±110 


150 


4843±35 


470±35 


1505±70 


6242±110 


100 


4802±35 


511±35 


1464±70 


6312±110 



Table II. 



AgcD(MeV) 


a(m b ) 


<B 2 > 


fH 


m cr (GeV) 


(/^^(MeV) 


200 


-0.19 


0.10 


-0.11 




260 


150 


-0.17 


0.19 


-0.08 


90 


210 


100 


-0.15 


0.41 


-0.12 


160 


210 



Table III. 



ÅQCD(MeV) 


l/a cc , (MeV) 


l/a bc , (MeV) 


l/a 66 , (MeV) 


200 


630 


790 


1240 


150 


540 


700 


1120 


100 


450 


590 


980 



Table IV. 



A QCD (MeV) 


a(2/j bc ) 


<B 2 > 


fH 


/s c (MeV) 


200 


-0.24 


0.35 


-0.44 


370 


150 


-0.22 


0.74 


-0.34 


540 


100 


-0.19 


1.93 


-0.54 


780 



Table V. 



ÅQCD(MeV) 


200 


150 


100 


B c -T 


-0.10 


-0.08 


-0.14 



37 



